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Temporal graphs are useful tools to model dynamic network topologies found in many 
applications. In this paper, we address the problem of constructing a convergecast tree 
on temporal graphs for data collection in dynamic sensor networks. Two types of con- 
vergecast trees, bounded arrival time convergecast tree and minimum total arrival time 
convergecast tree are defined as useful structures for low delay data collection. An 
O(T(n+m)) time centralized algorithm is proposed to construct a bounded arrival time 
convergecast tree, where n is the number of nodes, m is the number of edges, and T is the 
lifetime of the given temporal graph. The algorithm presented is an offline algorithm and 
assumes that all information about change in the graph topology is known apriori. It is 
then shown that the problem of constructing a minimum total arrival time convergecast 
tree is NP-complete, and an O(T (n + m)) time centralized, offline heuristic algorithm 
is proposed to address it. The heuristic algorithm is evaluated with experiments on four 
real life data sets. 


Keywords: Convergecast tree; temporal graph. 


1. Introduction 


Temporal graphs [20] are useful tools to model dynamic network topologies where 
the edge set and/or node set vary with time. Examples of such dynamic topologies 
can be found in many practical scenarios such as in delay tolerant networks (DTN) 
[19], vehicular ad-hoc networks [II], mobile ad-hoc networks (MANET) [13] [14] 
and social networks [28]. Several works on sensor networks have exhibited dynamic 
network topologies such as vehicular sensor networks to monitor urban areas [21] 
[23], DTN based systems for agricultural data collection [25], Bikenet [10] etc. In 
such systems, an end-to-end path between two nodes may not exist at any time. 
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Convergecast trees have been used to address the data aggregation problem in 
static sensor networks to collect data from the sensors (nodes of the tree) to a base 
station (the root node) [i], [29]. In a static network where the topology does not 
change, a spanning tree can be used both as a convergecast tree and a broadcast 
tree. However, this is not true for a temporal graph. In a temporal graph, an edge 
can only be used for traversal if it exists at the time of the traversal. A path from 
a node u to a node v is said to exist if there exists a sequence of edges from u to 
v that exist in strictly increasing order of time. Such paths have been referred to 
in the literature as time respecting path or journey [6]; we will use the term 
journey in the rest of this paper for such paths. It is obvious that the existence of a 
journey from u to v in a temporal graph does not necessarily imply the existence of 
a journey from v to u. Thus, a broadcast tree and convergecast tree may be different 
for temporal graphs; it is even possible that one exists while the other does not. 
Several recent works have addressed the problem of constructing broadcast trees 
[17], [I8]. However, there does not exist any work on constructing convergecast trees 
in a temporal graph. 

A journey from a node u to v in a temporal graph is said to be the foremost 
journey if it has the least arrival time at v among all such journeys from u to v. In 
this paper, we first define two types of rooted convergecast trees for temporal graphs, 
Bounded Arrival Time Convergecast Tree (BATCT), and Minimum Total Arrival 
Time Convergecast Tree (MTATCT). In a bounded arrival time convergecast tree, 
the maximum arrival time of a journey from any non-root node to the root node is 
bounded by the maximum of arrival times of the foremost journey from any non- 
root node to the root node in the given temporal graph. Similarly, in a minimum 
total arrival time convergecast tree, the sum of the arrival times of all journeys from 
all non-root nodes to the root node is the minimum among all such convergecast 
trees. Such convergecast trees can be useful for data collection with low delays 
in dynamic sensor networks. We present a O(T (n + m)) time centralized, offline 
algorithm to construct a BATCT, where n is the number of nodes, m is the number 
of edges in the temporal graph, and 7 is the total time for which the temporal graph 
exists. We next show that the problem of computing a MTATCT is NP-Complete, 
and present a O(T (n + m)) time centralized, offline heuristic algorithm for it. The 
heuristic algorithm is also evaluated on four real life data sets and it is shown that 
the algorithm performs reasonably well for most of the cases. To the best of our 
knowledge, no algorithm to construct any type of foremost convergecast trees has 
been proposed before this work. 

The rest of this paper is organized as follows. Section [2] discusses some relevant 
related work in the area. SectionB]describes the details of the model for the temporal 
graph and other assumptions made to address the problem. Section [4| formally 
defines the problem and terminology related to temporal graphs. Section[5] presents 
a centralized algorithm to compute foremost journeys from all non-root nodes to 
the root node. Section [6] presents a centralized offline algorithm for the BATCT 
problem. Section[7]proves that the MTATCT problem is NP-complete and describes 
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a centralized, offline heuristic to solve the problem. It also provides the details of 
experiments done to evaluate the proposed heuristic. 


2. Related Work 


Problems related to paths and spanning trees in temporal graphs have been 
addressed by researchers in recent times. In a temporal graph, a path from a node 
u to anode v is a sequence of edges from u to v such that the edges in the sequence 
exist in strictly increasing order of time. Paths in a temporal graph have been 
referred to as time respecting path [24], journey [6] etc. In Ref. [4], Bui-Xuan et al. 
have defined and proposed algorithms to find three different types of journeys in 
temporal graphs, namely foremost journey, fastest journey, and shortest journey. 
Wu et al. [B0] have defined another type of journey called latest departure time 
journey. They have proposed algorithms to compute foremost, fastest, and shortest 
journey from a node to all other nodes, and latest departure time journey from all 
nodes to a given node. Some other works that address path related problems in 
temporal graphs are [9], [26], and [BI]. 

In works related to broadcast trees, Huang et al. [I7] have defined two types of 
broadcast trees for weighted temporal graphs, minimum spanning tree with earliest 
arrival times and minimum spanning tree with smallest total weight. A polynomial 
time algorithm is proposed for the first problem. The second problem is shown to 
be NP-complete and an approximate solution is proposed. Ikuta et al. [18] have 
proposed an integer linear program based solution for the minimum spanning tree 
with smallest total weight problem and experimentally shown that their approach 
performs better than the solution proposed in Ref. [I7]. 

As discussed earlier, the existence of a broadcast tree does not imply the exis- 
tence of a convergecast tree in a temporal graph. Recently the problem of data 
aggregation in temporal graphs has received considerable attention where converge- 
cast trees can be useful. In Ref. [8], Cornejo et al. have addressed the problem 
of accumulating tokens distributed among the nodes in a distributed environment 
to a single node in the system. It is shown that no data aggregation algorithm 
can be competitive with offline algorithms for every dynamic graph, and an algo- 
rithm which guarantees good competitive ratio when the graph is well connected 
is proposed. Bramas et al. have proposed an online algorithm to address the 
data aggregation problem for dynamic graphs in the presence of oblivious, adap- 
tive, and random adversaries in distributed environments. In Ref. [3], Bramas et al. 
have proposed a foremost convergecast tree based approach to solve the minimum 
data aggregation schedule problem for dynamic wireless sensor networks where they 
have represented the dynamic graph as evolving graph. They have assumed that a 
node can receive data from only one neighbour at each timestep; our model allows 
multiple messages to be received at a node in a single timestep. To compute the 
minimum data aggregation schedule, the arrival time of the foremost journey with 
the maximum arrival time in the foremost convergecast tree is computed. However, 
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no algorithm to construct the foremost convergecast tree for a given temporal graph 
is presented. 

In this paper, we address the problems of construction of a BATCT and a 
MTATCT. To the best of our knowledge there is no work in existing literature which 
proposes algorithms to construct such types of convergecast trees. 


3. System Model 


We represent a temporal graph by the evolving graphs [I2] model. In this model, 
a temporal graph is represented as a finite sequence of static graphs, each static 
graph being an undirected graph representing the graph at a discrete timestep. The 
total number of timesteps is called the lifetime of the temporal graph. A temporal 
graph is denoted by G(V,€) = {Go(V, Eo), Gi(V, €1),..., @r_-1(V, ET-1)} where V 
is the node set, E = Ue Ei is the edge set, and F is lifetime of the temporal 
graph. Each G; denotes the graph at timestep 7 with node set V and edge set 
Ei. The temporal graph G is available for time interval [0, 7). Temporal graph G 
can be equivalently represented by specifying the time intervals for which each 
edge is available. Thus an edge e € E between nodes u and v can be represented 
as e(u,v, (51, fi), (S2, f2),---, (Sk, fk)), where u,v € V,u Æ v, and a pair (s;, fi), 
1<i<k indicates that the edge exists for the time interval [s;, fi), where 0 < s; < 
fi < T (and hence the edge exists in all the static graphs G,,Gs,41,.-.,G@y,-1). 
Also, if an edge e has two such pairs (s;, fi) and (sj, fj), si A sj and if si < sj 
then f; < sj. An edge at a single timestep is called an instance of that edge. For 
simplicity, in the rest of this paper, when the exact time intervals for which an edge 
e exists are not important, we denote the edge e between nodes u,v E€ V by ew. An 
instance of ew at time t is denoted as et. Also, we use n and m to denote |V| and 
|E| respectively. 

For this work, we have assumed that the node set remains unchanged throughout 
the lifetime of the temporal graph, and only the edge set changes with time. Also, 
for simplicity it is assumed that all information about the changes of the edge set is 
known a priori. Note that while this assumption may not be valid for non-periodic 
temporal graphs with lifetime 7, these algorithms may still be applicable to a single 
period (of duration 7) of a periodic graph after learning the edge change pattern 
from past periods. The traversal time of each edge is taken to be one timestep. 
During traversal, waiting at a node is allowed, i.e., a traversal following a journey, 
on arrival at a node u, can wait at u for any duration of time before traversing the 
next edge in the journey. We also assume that there is at least one journey from all 
non-root nodes to the root node. This assumption ensures the existence of at least 
one convergecast tree in the given temporal graph. 


4. Problem Definition 


In this section we formally define the problems addressed in this paper. The input 
for the problems is a temporal graph G(V,€) and a predefined root node, r € V. 
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The lifetime of G is T. The underlying graph of G is defined as the graph Gu (V, Eo U 
E1 U... UEr—1); thus, it is the union of the static graphs at all timesteps. We first 
define the following terms which will be used to define the problem. 


Definition 1. Journey [6]: For a given temporal graph G(V,€), a journey 
Jtte*1 (u,v) starting at timestep t from a source node u € V and arriving at a des- 
tination node v € V at timestep t + 1, is a sequence {(e1, t1), (e2, t2),.-., (ek, th) } 
where e; € E, and (e1, €2,..., ex) is a walk in Gy with u as the starting node of e1 
and v as the ending node of ex, tı > t, ti < T, e; exists at t;, and (ti + 1) < ti41. 
Timestep t is called the departure time and timestep tp +1 is called the arrival time 
of the journey. 


Note that since waiting at a node is allowed, the first edge traversed by a journey 
from u with departure time t can be an edge instance existing at a timestep tı 
greater than t also (i.e. after waiting at u for the duration tı — t). 


Definition 2. Length/Hop count of a journey [4]: The length or hop count of 
a journey J‘s**# (u,v) is the number of edges in the journey. 


Definition 3. Foremost journey at timestep t [4]: A foremost journey from a 
source node u € Y to a destination node v € VY with departure time t is the journey 
with the earliest arrival time at v among all journeys from u to v with departure 
time t. 


It may be noted that if Jets (u,v) and jee (u,v) are foremost journeys from 
a node u to a node v with departure times tı and tə respectively, then if tə > tı, 
then tr, < tf (as Jı can wait at u for (t2 — tı) time and traverse along J2 to arrive 
at ty, at v). In the rest of this paper, if t is not explicitly specified, it refers to the 
foremost journey from u to v with the earliest arrival time at v among all foremost 
journeys from u to v with any departure time. Also note that as waiting at any node 
is allowed, a foremost journey with departure time 0 is always a foremost journey. 


Definition 4. Latest departure time [4]: Latest departure time from a source 
node u € V to a destination node v € VY is the latest timestep when a journey can 
start from u and reach v within the lifetime of the temporal graph. 


Definition 5. Journey induced subgraph [5]: Let there be a temporal graph 
G(V,€) with lifetime T and a set of journeys J = {J1, J2,...,In}, with V(J;) as 
the node set and €(.J;) as the edge set of each J;. Then, G7(U; V(Ji), U; E(J;)) is a 
journey induced subgraph of G. 


Definition 6. Convergecast tree: For a given temporal graph G(V, E£) and a node 
r € V, if T is a journey induced subgraph of G such that there exists exactly one 
journey from all x € V — {r} to r then T is a convergecast tree rooted at r for G. 


In a convergecast tree each node has a set of parent nodes and a child node 
(except the root node). For a node v, the set of parent nodes is denoted by P(v); 


390 S. Mandal & A. Gupta 


Fig. 1. (a) A temporal graph, (b) bounded arrival time convergecast tree rooted at node r, 
(c) minimum total arrival time convergecast tree rooted at node r. 


note that for leaf nodes P(v) = Ø. The child node of v is denoted by c,; for root 
node r, cr = ġ. Note that even though the existence of a convergecast tree does 
not imply the existence of a broadcast tree, a broadcast tree in a temporal graph 
is the same as a convergecast tree with each edge instance et, of convergecast tree 
replaced with the edge instance ene forO0<t<7-1. 


v 


Definition 7. Bounded arrival time convergecast tree (BATCT): Given a 
temporal graph G(V, E) and a node r € V, a convergecast tree T is a bounded arrival 
time convergecast tree rooted at r if the maximum of arrival time of any foremost 
journey from any node x € V — {r} to r in T is less than or equal to the maximum 
of the arrival times of foremost journeys from any node x to r in G. 


Definition 8. Minimum total arrival time convergecast tree (MTATCT): 
Given a temporal graph G(V,€) and a node r € V, a convergecast tree T is a 
minimum total arrival time convergecast tree rooted at r if the sum of the arrival 
times of the journeys in T from all nodes « € V — {r} to r is minimum among all 
such convergecast trees. 


For the temporal graph shown in Fig. Ma), J*®(f,r) = {((f,9) ,4), ((g,r),5)} 
is a journey from node f to node r with departure time 4, arrival time 6 and 
hop count 2. This is also the foremost journey from node f to node r. The lat- 
est departure time for journeys from node f to node r is 6. A BATCT rooted 
at node r is the structure shown in Fig. [I{b) for which the journeys with the 
maximum arrival time are J°°(e,r) = {((e,c), 0), ((c, b), 6), ((b,r),8)}, J19 (d, r) = 
{((d, c), 1), ((c, b), 6), ((,7),8)}, J°P(e,r) = {((c,b),6), ((0,r),8)}, J59 @,r) = 
{((b,r),8)} with arrival time 9. In the original graph the maximum arrival time of 
any foremost journey to node r is J”? (e,r) = {((e, c), 0), ((c, b), 6), ((b, 7), 8) } with 
arrival time 9. Figure [{c) shows a MTATCT for the given temporal graph with 
total arrival time 49. 

In this paper, we address the problems of constructing bounded arrival time 
convergecast tree and minimum total arrival time convergecast tree for a given 
temporal graph with a given root. To solve these problems, we use an algorithm 
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that, for each node v and each timestep t, computes the foremost journey J from v 
to the root node r with departure time t if the arrival time of J at r is smaller than 
the arrival time of the foremost journey from v to r with departure time tı for all 
tı >t. This algorithm is described in the next section. 


5. Finding All Foremost Journeys 


For any node v, the foremost journeys from v to the root node r with departure 
time t for all t < T can be found in O(nm7T (log(7) + log(n))) time using the 
algorithm described in Ref. [4]. We propose a faster algorithm which computes 
foremost journeys from all non-root nodes to the root node r in O(T (n+ m)) time. 
To compute this, we process each edge instance of G in non-ascending order of their 
time of existence. If edge instances of multiple edges exist at the same timestep, 
they can be processed in any arbitrary order. Consider a non-root node v and a 
timestep t. The proposed algorithm computes a foremost journey J from v to r 
with departure time t, if the arrival time of J is smaller than the arrival time of the 
foremost journey from v to r with departure time tı for all tı > t. 

Existence of an edge instance et „ denotes the existence of the journey J¢‘+1(u, v) 
and J‘é‘+!(y,u). If there is an edge instance et, and a journey J" 'f(u,w) (or/and 
J‘utt(v,w)) such that tı > t then a journey Jt" (v, w) (or/and J's (u, w)) also 
exists. Following this idea, we compute journeys with departure time ¢ from all 
non-root nodes to the root node by using journeys with departure time tı > t. We 
process edge instances in non-ascending order of their time of existence starting with 
timestep T — 1. It terminates when all edge instances of 0t” timestep are processed. 

The pseudocode is given in Algorithm [I] The algorithm takes a temporal graph 
G(V,€) with the root node r € VY as input. A[v][ét], N[v][t], and Hf[v][t] are 2- 
dimensional arrays maintained to store the arrival time, the next node of v, and 
the hop count of the foremost journey from node v to r with departure time t 
respectively. Afv][t] is initialized to oo for all v, t. Note that A[v][t] is updated only 
if the arrival time of the foremost journey from v to r with departure time t is 
smaller than the arrival time of the foremost journey from v to r with departure 
time tı for all t > tı; otherwise A[v][t] remains as oo and the values of N[v]|t] and 
H{v|[t] are not relevant. We also maintain an array F'[v] to store the departure time 
of a foremost journey from node v to r (v Æ r) computed till the current timestep. 
Another array PF|v] stores the departure time of the foremost journey from v to r 
with the minimum arrival time at r among all foremost journeys from v to r with 
departure time t for all t > F'[v]. If more than one foremost journeys with departure 
time t > Fv] have the same arrival time at r, PF'[v] stores the one with the latest 
departure time. 

Algorithm [I] shows the pseudocode for processing the edge instances in non- 
ascending order of timesteps of their existence. While processing an edge instance 
et, Algorithm[]in Lines [4] or[6] first checks if either u or v is the root node r. If u is 


Uv? 


the root node, then the algorithm updates arrival time (A|v][t]), hop count (H[v][#]) 
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Algorithm 1 computeAllJourneys(G(V, E), r) 
Input: A temporal graph G(V,€), with root node r, lifetime T. 
Output: Foremost journey from each node to r for each departure time 
1: Initialize: Vt; < T N[r][ti] := r, Alr][ti] := ti, H[r][ti] := 0 
2: Initialize: Vv; E€ V — {r}, Vti < TN [vi] [ti] := 6, Alvi] [ti] := 0, H[vi][ti] := co 
F[y;] := 00, PF [vi] := 00 
3: for all ec’, t = (T — 1) +0 do 


Uv? 


4 if v=r then 

5 update Value(u, v, t+1, t) 

6: else if u = r then 

7 update Value(v, u, t+1, t) 

8 else 

9: if F'[u] = œ and Fv] Æ œ then 

10: if Flv] At then 

11: update Value(u, v, Ffu], t) 

12; else if PF|v] 4 œ then 

13: update Value(u, v, PF/v/, t) 

14: end if 

15: else if Flu] 4 œ and Fv] = œ then 
16: if Flu] At then 

17: update Value(v, u, Flu], t) 

18: else if PF|u] 4 oo then 

19: update Value(v, u, PFfuj], t) 
20: end if 
21: else 
22: if Alu][F[u]] > A[v][F[v]] then 
23: if F[v] >t+1 then 
24: update Value(u, v, Ffu], t) 
25: else if PF|v] 4 œ and Aļu]|F[u]] > Al[v][PF[v]] then 
26: update Value(u, v, PF[v], t) 
27: end if 
28: else if A[v][F[v]] > Aļu][F[u]] then 
29: if F[u] >t+1 then 
30: update Value(v, u, Ffuj, t) 
31: else if PF|u] 4 co and Aļ[v]|F|v]] > Alu][PF [uj] then 
32: update Value(v, u, PFfu], t) 
33: end if 
34: end if 
35: end if 


36: end if 
37: end for 
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and next node (N [|v][t]) of the foremost journey from v to r with departure time t, 
with t+ 1, 1 and u respectively. Similar action is taken if v is equal to r. If none of u 
or v is the root node r, then the algorithm checks if only one of u or v has a journey 
to r (Lines Dļland [5). If only u has a journey to r, then the algorithm updates the 
arrival time (A[v][t]) of the foremost journey from v to r with departure time t, 
with the arrival time of the foremost journey from u to r with departure time Ffu] 
(if Flu] > t+ 1 or with departure time PF [u] unless PF'[u] = 00). It also updates 
the hop count (H[v][t]) of the foremost journey from v to r with departure time t, 
to one plus the hop count of the foremost journey from u to r with departure time 
Flu] (if Flu] > t+ 1 or with departure time PF Ju] unless PF[u] = oo) and sets u 
as the next node (N[v][t]) to v. Similar action is taken if only v has a journey to r. 
If none of these conditions are satisfied, then both u and v have journeys to r. In 
this scenario, Algorithm [I] checks if the arrival time (A[u][F[u]] if F[u] > t+ 1, else 
Alu][PF [u]] unless PF[u] = 00) of any foremost journey from u to r with departure 
time tı for all tı > t is smaller than the arrival time of the foremost journey from 
v to r (Lines [22] and 28). If this condition is true, then the algorithm updates the 
arrival time (A[v][t]) of the foremost journey from v to r with departure time t, 


with the arrival time of the foremost journey from u to r with departure time F'[u] 
(if Flu] > t+ 1 or with departure time PF[u] unless PF'[u] = oo). It also updates 
the hop count (H[v][t]) of the foremost journey from v to r with departure time t, 
to one plus the hop count of the foremost journey from u to r with departure time 
Fu] Gf Flu] > t+ 1 or with departure time PF[u] unless PF'[u] = oo) and sets 
u as next node (N{v][t]) to v. Similar action is taken if the arrival time (A[v][F'[v]] 
if Flv] > t+ 1, else A[v][PF[v]] unless PF[v] = co) of the foremost journey from 
v to r with departure time tı for all tı > t is smaller than the arrival time of the 
foremost journey from u to r. It terminates when all edge instances of timestep 0 
are processed. 

As an example consider the temporal graph shown in Fig. [I{a), with root node 
r. The algorithm first processes edge instances e?,, e®,. It adds the foremost journey 
from a with departure time 9 and discards e®,, as it does not lead to any journey 
to r. Then similarly it processes e., e8,, e8, and discards e°,. Edge instance e8, 
satisfies the condition stated at Line] or [5] and the algorithm adds the foremost 
journey from c via a with departure time 8. Again êen satisfies the condition stated 
at Line P2]or B8]and as b has a journey to r with lower arrival time, the algorithm 
adds foremost journey from node c. It terminates when processing of e}, is done. 


Theorem 9. Let G(V,E) be a given temporal graph with the root node r and life- 
time T. For each non-root node v € V and each timestep t, Algorithm U] correctly 
computes the foremost journey J®* (v, r) from v to the root node r € V with depar- 
ture time t if the arrival time of J**t (v,r) atr is smaller than the arrival time of 
the foremost journey from v to r with departure time tı for all tı >t. 


Proof. In Algorithm [I] edge instances are processed in non-ascending order of 
timesteps of their time of existence. We have to show that for all non-root nodes 
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Procedure 2 updateValue(target, source, stime, time) 
N{target]|time] := source 
Altarget][time] := Alsource]|stime] 
H [target] [time] := H{source][stime] + 1 
if F'[target] > time then 
PF target] := F [target] 
end if 
F [target] := time 


v E€ V— {r} at any timestep 0 < t < T, if there exists a foremost journey J‘? (v, r) 
from node v to r with arrival time smaller than the arrival time of the foremost 
journey from v to r with departure time tı for all tı > t then Algorithm [I] correctly 
computes it. We prove this by induction on departure time t of foremost journeys. 

Base Case: t = T — 1: Let v be a neighbour of r and let there be an edge 
instance e771. Then Line B] or [7] of Algorithm [I] inserts it as a foremost journey 
from v to r with departure time 7 — 1. In this scenario the condition in Line [4] or 
Line [6] gets satisfied and the values of Alv][t] := t + 1, N[v|[t] := r, H[v][t] := 1 and 
F |v] := t get updated by Procedure B] As t = T — 1 the value of PF[v] remains oo. 
As traversal time of each edge is 1 timestep, the existence of another journey with 
earlier arrival time with departure time 7 — 1 is not possible. Hence this shows that 
Algorithm [I] correctly computes the foremost journeys with departure time 7 — 1. 

Inductive step: We assume that Algorithm [J] correctly computes the foremost 
journeys with departure time t > (7 —k), ts := (J —k). We have to show that Algo- 
rithm GJ correctly computes the foremost journeys with departure time ts — 1. Since 
edge instances are processed in non-ascending order of their time of existence, then 
before processing edge instances of timestep ts — 1, all edge instances of timesteps 
ts to T — 1 are processed. Before computing the foremost journeys with departure 
time ts — 1 we have all the foremost journeys with departure times t, to T — 1. If 
any node v, has a journey J* 71t (v,r), then J needs to go through some other 
journey J¢s:¢7 (u,r) such that t, > ts. Using this property Algorithm [I] computes 
the possible journeys from node v to r with departure time ts — 1. Let efs! be 
the edge instance such that, there is a journey with departure time Ffu] > ts (or 
PF \u| > ts if Flu] = ts — 1 and PF[u] 4 oo) from u to r with earlier arrival time 
than the arrival time of the journey with departure time Fv] > ts — 1 from v to 
r. Then the condition at Line or is satisfied. The algorithm computes the 
foremost journey from v to r with departure time ts — 1 which is considering ets ', 
by comparing A[v][F[vu]] with A[u][F [ul] if Flu] > ts, or A[v][F[e]] with Afu][PF [ul] 
if F[u] = ts and PF[u] # co. As the condition stated in Line [f] in Procedure B] 
ensures that either only PF Ju] or both F'[u] and PF[u] store timestep greater than 
or equal to t, and all information stored in N, A, H, F and PF for timestep t, to 
T — 1 are correct then the computed foremost journey with departure time ts — 1 
from node v to r is correct. This proves the inductive step. oO 
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Theorem 10. The time complexity of Algorithm] is O(T (n + m)). 


Proof. In Algorithm[]] the initialization phase initializes the 2-dimensional arrays 
A, N, and H and the 1-dimensional arrays F and PF. The arrays A, N, and H each 
has nT entries. The arrays F and PF each has n entries. Thus the initialization 
step takes O(n7) time. At each step the algorithm processes an edge instance, and 
this processing of a single edge instance takes O(1) time. The maximum number of 
edge instances is m7. Hence the running time of the algorithm is O(7(n + m)). O 


6. Constructing a Bounded Arrival Time Convergecast Tree 


In this section, we present a centralized algorithm to solve the problem of construct- 
ing a BATCT for a given temporal graph G(V,€) and a given root node r € V. 
The algorithm to construct the BATCT works in two phases. In the first phase, 
we find the foremost journey J from each non-root node v to the root node r with 
each departure time t, if the arrival time of J at r is earlier than the arrival time 
of the foremost journey from v to r with departure time tı for all tı > t. This is 
done using Algorithm[] given in Sec. [5] The second phase starts with an empty con- 
vergecast tree T and iteratively adds nodes to T. At each iteration, the algorithm 
considers the set of nodes that have not yet been added to the convergecast tree. 
From the set of foremost journeys of all nodes in this set of nodes, the foremost 
journey with the maximum arrival time is added to T. Note that a journey can add 
multiple new nodes to T. The process is repeated until all nodes in V are added to 
T. The pseudocode of the algorithm is shown in Algorithm[3] Various symbols used 
in Algorithm B] are listed below. 


(1) maaForemost(S): returns the node which has the foremost journey with the 
maximum arrival time from a set of nodes S. 

(2) startForemost(v): returns the timestep at which the foremost journey from node 
v to r traverses the edge instance incident on v. 

(3) nextNode(J':'f(s,d),u,ti): returns the next node after u in the journey 
J'=*f (s d), where u is an intermediate node in J‘:'f(s,d) that leaves u at 
timestep tı such that tı < tr. 

(4) v.visited: set to 1 when a node v is added to T; initially it is set to 0. 

(5) nextTime(J':"'f (s, d), w,t’): returns the timestep when the journey Jtet (s, d) 
departs from an intermediate node w after arriving at w at timestep t. 


As an illustration, consider the temporal graph shown in Fig. [{a) with node 
r as the root node. The foremost journeys from nodes c,d,e all have arrival 
time 9, the foremost journey from b has arrival time 4, the foremost journeys 
from a,g have arrival time 3, and the foremost journey from f has arrival time 
6. The algorithm randomly chooses a node from c,d,e and adds the foremost 
journey from the node to the tree. Suppose this node be d and the journey 
J=? (d, r) = {((d,c), 1), ((c, b), 6), ((b, 7), 8)} is added to the tree. In the next step, 
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Algorithm 3 boundedTimeConvergecastTree(G(V, E), r) 
Input: A temporal graph G with node set V edge set E, lifetime 7 and root r. 
Output: A bounded arrival time convergecast tree T(V’, E’), V’ = V rooted on r. 
1: computeAllJourneys(G(V, E), r) 
2: Initialize: S := V — {r}, V' := {r}, r.visited:= 1, P(r) :=0, cr := o 
3: while S 4 ý do 
4: v := mazForemost( S) 
5 t := startForemost(v) 
6: J‘*'t(v,r) := the foremost journey from v to r at time t 
T: while v 4 r do 
8 v.visited := 1 
9: w := neatNode( J's (vu, r), v, €) 
10: P(w) := P(w) Uv 


11: Cy (= Ww 

12: V :=V'U {v} 

13: E' := E' U {e(v, w, (t,t + 1))} 
14: S:=S\v 

15: if w.visited = 0 then 

16: v := w 

17: t := nextTime(J'"! (v,r), w,t + 1) 
18: else 

19: break 

20: end if 

21: end while 


22: end while 
23: return(T(V’, E’)) 


it adds J” (e,r) = {((e,c), 0), ((c, b), 6), ((b, 7), 8)}. At this step, after adding edge 
instance e? , it reaches a node c for which c.visited = 1 and it stops there 
and selects the next journey which is J*®(f,r) = {((f,g),4), ((g,7),5)}. Finally 
J?3(a,r) = {((a,r),2)} is added. The resultant tree is shown in Fig. [{b). 


Lemma 11. Let A a (v,w) and je (u,w) be two foremost journeys in a tem- 
poral graph G(V, E), such that ty, < tfa. Letx € V be the first node which is common 


in Jı and J2. Then there exists a journey Jh (v, w). 


Proof. Let Jp" (x, w) and Jt" (x, w) be the journeys followed by Jı and J2 
respectively from x to w. We first prove that Jı reaches x earlier than J2. There 
are two possible cases. 


(1) tp, < t2: In this case, it is obvious that Jı reaches x earlier than J2. 
(2) tp, > t2: We prove this by contradiction. Let Jz reach x earlier than Jı. This 
implies that there is a journey from u which follows Jp (u, w) upto node z, 
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then follows J 1 (x, w) and reaches w at ty,. It is given that tf < tf. This 
contradicts the fact that Jz is the foremost journey from u to w. Hence Jı 
reaches x earlier than J2. 


As Jı reaches x earlier than J2, there exists another journey ioe (v,w) which 
can reach w by following JE stay 
reach w. o 


v, w) upto node z, then following J; z, w) to 


Lemma 12. Let ae (v,r) and i (u,r) be two foremost journeys from two 


non-root nodes v and u to the root node r with departure times tı and tz respectively 
computed by Algorithm [E] Let x E€ V—r be the first node which is common in Jı, 
J2 and Jı, Jz reach at x at t3, t4 respectively where t4 > t3. Then both Ji, Jo follow 
same journey I (x,r) such that ts > t4. 


Proof. We prove this lemma by contradiction. Let Jı, Jo follow two different jour- 
neys sists (x,r) and eats (x,r) respectively. It is clear from the definition of a 
journey that t4 > t4 and ts > ts. As both J4 and Js arrive at r at ty, then J4 # Js 
implies that for node x, Algorithm[]stores information about two different journeys 
from x to r with the same arrival time ty. This contradicts the fact that, while com- 
puting the foremost journey from any non-root node zx to r, if two or more foremost 
journeys with different departure times have the same arrival time ts, Algorithm [] 
stores the foremost journey with the latest departure time. Hence J4 = Js and the 
journey from x to r must have departure time ts > t4 (otherwise Jz (if t4 > ts > ts) 
or both Jy, J2 (if t3 > ts) is invalid). O 


Theorem 13. Algorithm [J] correctly computes a BATCT, T rooted at node r in 
O(T(n + m)) time. 


Proof. Algorithm [B] computes the foremost journeys from all non-root nodes to 
the root node r using Algorithm [] This takes O(7 (n + m)) time. 

At each iteration, Lines [4] to of the algorithm adds the foremost journey 
Jt=ttf (v, r) such that v is not in T and the arrival time of J‘='f (v, r) is the maximum 
among foremost journeys from all remaining nodes. While adding a journey, after 
addition of each node, Line [4] removes the added node from the list of remaining 
nodes. The while-loop from Lines [7] to PI] stops if next node is r or a node w such 
that w € T. Line [5] checks whether w is already present in T or not. This ensures 
that each node is added only once to T. As w is already present in T then there 


is a journey gents (w,r), tp < t in T. Lemmas [I] and [2] ensures that there is 
a journey A (v,r) via w. This is true for any node v. This shows that in T, 
there is a journey from each node in VY — {r} to r with arrival time bounded by the 
maximum of the arrival times of foremost journeys from any non-root node to r in 
the temporal graph. This procedure also ensures that in T each node has exactly 
one journey to r. 
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Sorting the nodes according to the arrival times of their foremost journeys takes 
O(n + T) time as arrival time of any journey is an integer and bounded by T. 
Each node addition takes O(1) time and each node is added only once. Hence this 
process of adding the nodes into T takes a total of O(n) time. So running time of 
the algorithm is O(T (n + m)). Oo 


7. Constructing a Minimum Total Arrival Time Convergecast Tree 


We first prove that construction of MTATCT for a given temporal graph is a NP- 
Complete problem. 


Theorem 14. The problem of constructing a minimum total arrival time converge- 
cast tree is NP-Complete. 


Proof. We first define a decision version of the problem, referred to as D-MTATCT, 
as follows: Given a temporal graph G(V, E), a root node r € V and a positive integer 
j, does there exist a convergecast tree T(V,€’) such that €’ C E and the sum of 
arrival times of all journeys from nodes in (V — {r}) to r is at most j? 

We first show that D-MTATCT is in NP. Consider a certificate ((G(V, E€), r, j}, T) 
where G(V,€) is a given temporal graph, r € V is the root node, j is the required 
upper bound of the sum of arrival times of all journeys from nodes in (V — {r }) to r, 
and T is the given MTATCT on G. Given the certificate, we can check whether the 
sum of arrival times of all journeys in T from each node in (V — r) to r is at most j 
in polynomial time. To perform this check, the arrival time of the journey from each 
node to r is computed. This is done by starting from each node and following the 
available edges to the child node of the current node at each step until the journey 
arrives at r. Whether E’ C E, can also be checked easily in polynomial time. Hence, 
D-MTATCT is in NP. 

Next we prove that there exists a polynomial time reduction from the set cover 
problem [I5], which is known to be NP-complete, to the MTATCT problem. 

Consider an instance of the set cover problem (U,S,k), where U = 
{a1,d2,...,@n} is a set of n elements, S = {S1,S2,...,Sm}, a collection of m sets 
such that Us,es Si =U, and a given integer k. The set cover problem is to decide 
if there exists a subset, S’ C S of S such that Uses S; =U and |S"| < k. 

From the given instance of the set cover problem, we construct an instance of 
D-MTATCT problem. We construct a temporal graph G(V, E) as follows: 


e A node S; is added to VY for each set S; € S. 
e A node a; is added to V for each element a; € U. 
e A root node r is added to V. 

e An edge e(a;, Si, (0, T)) is added to £ if a; € Sj. 
e VS; an edge e(r, Si, (0, T)) is added to £. 
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Fig. 2. Temporal graph construction from set cover problem instance. 


Hence the temporal graph G has n +m + 1 nodes. The traversal time of each edge 
is set to 1 timestep, j = 2(k +n) + (m — k) and lifetime T = 2. An example of the 
temporal graph construction is shown in Fig. B] 

We first show that, if there is a solution for the D-MTATCT problem instance 
then we have a solution for the set cover problem instance. For a given convergecast 
tree T rooted at node r for G, such that the sum of arrival times of journeys from 
all nodes in (V — {r}) to r is 2(k +n) + (m — k), we construct a set cover S C S 
such that |S’| = k. Consider the set S” = {S;| node S; has journey to r with arrival 
time 2 in T}. We show that S’ covers all elements in U and |S’| = k. We prove 
this by contradiction. Let us assume that a; is an element such that a; E€ U and 
aiU sies Si- AS S’ includes all S; whose corresponding S; nodes have journey to 
r with arrival time 2 then our assumption can be true only if one of the following 
conditions hold. 


(1) Node corresponding to a; is not included in T. This contradicts the fact that T 
is a convergecast tree on G. 

(2) There is a journey J from node corresponding to a; to r with arrival time 1. As 
no edge between node a; and r is in €, existence of J with arrival time 1 implies 
that J arrives at r at timestep 1 from a; via some S; node. This contradicts the 
fact that traversal time of each edge is 1 timestep. 


As both the conditions lead to contradiction, Us. ¢5, Si = U. To prove |S"| = k, two 
cases are possible. If |S’| = 1 > k, this implies that, in T, there are l S; nodes which 
have journey to r with arrival time 2. As all a; nodes are at two hop distance from r 
in T and 7 = 2, the arrival time of each journey from each a; node to r is 2. Hence 
the sum of arrival times of all journeys in T is 2(1+n)+(m—l) > 2(k+n)+(m—k&), 
which is a contradiction. If |$’| = l < k, this implies that, in T, there are | S; 
nodes which have journey to r with arrival time 2. As 7 = 2, possible arrival times 
of each journey from each S; node are 1 and 2. The arrival time of each journey 
from each a; node to r is 2. Hence the sum of arrival times of all journeys in T 
is 20 +n)+(m-—l) < 2(k +n) + (m — k), which is again a contradiction. Hence 
|S’| =k. 
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Fig. 3. Bad case scenario for Algorithm [B] 


Next we show that, if there is no solution for the D-MTATCT problem instance 
then there is no solution for the set cover problem instance. To show this we show 
that if we have a solution for the set cover problem instance then we have a solution 
for the D-MTATCT problem instance. For a solution S’ C S, |S'| = k for the set 
cover problem, we can construct a convergecast tree T(V, E’) rooted at node r for 
G with total arrival time 2(k + n) + (m — k) in the following ways. 


(1) Add el s, to €' if a; € S; and S; € S$’. If a; is present in more than one set of 
S’ then we add edge instance with any one S; node. 

(2) Add els, to E’ if S; € S. 

(3) Add els, if S; ¢ 9”. 


We show that T is a convergecast tree for G with total arrival time 2(k + n) + 
(m — k). As S’ is a set cover, each node a; has at least one node S; € 9” as its 
neighbour. Only el. s; 1s added to T. For all S$; nodes, such that the corresponding 
set S; € S’, ebr is added to T. Hence all a; nodes have exactly one journey to r 
with arrival time 2 via some S; node such that a; € S; and S; € S’. Similarly all S; 
nodes corresponding to sets S; € S’ have only one journey to r with arrival time 2. 
All S; nodes corresponding to sets S; ¢ S’ have only one journey to r with arrival 
time 1. Hence total arrival time of the convergecast tree is 2(k +n) + (m-— k). O 


Algorithm B]can be used as an approximate solution for the MTATCT problem. 
However, in some cases, the ratio of the sum of arrival times of all journeys from all 
non-root nodes to the root node in the convergecast tree constructed by Algorithm [B] 
to the sum of arrival times of all journeys from non-root nodes to the root node 
in an optimum solution of the MTATCT problem may be high. For the temporal 
graph shown in Fig. [3| there are one hop journeys from all nodes uj, u2,...,Un—2 
to root node r with arrival time 1. There also exists the journey mee Cee r= 
{((Un—1, Un—2); 0), ((Un—2, un-3), 1), --., ((u1,r), T — 2)} where T = n, and the 
one hop journey JJ~'7 (un—1,r) = {((un-1;r), T — 1)}. Then the total arrival 
time of the MTATCT which includes Jz is n + T — 2; however the convergecast 
tree constructed by Algorithm [B] includes Jı and the total arrival time becomes 
(n — 1)(T — 1). Hence we use a different approach as described next. 

To address the MTATCT problem, we use Algorithm [I] that, for each node v 
and each timestep t, computes the foremost journey J from v to the root node r 
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Procedure 4 smallerHopCount(nodel, node2, time) 
1: if F[nodel] > time and H[node2][F'[node2]| > H[node1]|F|node1]] then 

2 updateValue(node2, nodel, F[nodel], time) 

3: else if PF'[nodel] 4 co and H|[node2]|F[node2]] > H[node1]|PF|node1]] then 
4 updateValue(node2, nodel, PF/nodel], time) 

5: end if 

6: if F[node2] > time and H[nodel]|F[node1]] > H|[node2]|F [node2]] then 

7: updateValue(nodel, node2, F[node2], time) 

8: else if PF'[node2] 4 co and H|nodel][F[nodel]] > H[node2]|PF'|node2]]| then 
9: updateValue(nodel, node2, PF [node2], time) 

10: end if 


with departure time t if the arrival time of J at r is smaller than the arrival 
time of the foremost journey from v to r with departure time tı for all tı > t. 
However, in this case, we modify Algorithm [[] such that if there are two foremost 
journeys from a node to the root node with the same arrival time then the algo- 


rithm stores the journey with the smaller hop count. To do this, while processing 


t 
Uv? 


an edge instance ef „, the algorithm checks if ef „ results in a journey from v and/or 
u to the root node which has arrival time equal to the arrival time of the cur- 
rent foremost journey from v and/or u. In this case the algorithm updates the 
journey information for v and/or u with the new journey if it has smaller hop 
count than the currently stored foremost journey. To perform this check we add 
another condition “else if A[v]|F[v]] = Alu][F[u]]” after Line B3)in Algorithm M] 
along with the existing conditions. The modified algorithm calls Procedure 4 if 
this condition is satisfied. We use this strategy because if a journey with smaller 
hop count from a node is added to the convergecast tree, then journeys with 
smaller arrival times can be chosen for larger number of remaining nodes. As an 
illustration consider the temporal graph G(V,€), shown in Fig. 4a). In G, there 
are two foremost journeys J;"(f,a) = {((f,e), 1), ((e, d), 2), ((d, c), 3), ((c, a), 4)}, 
J3?(f,a) = {((f,b), 0), (b,a), 4)} with the same arrival time 5 from f to a. While 
constructing a convergecast tree rooted on node a, if we store journey Jı as the 
foremost journey from f to a and construct a convergecast tree (shown in Fig.[4b)) 
then the total arrival time of all journeys from all non-root nodes to the root node 
is 21. However if we store Jo as the foremost journey from f to a, then this value 
for the convergecast tree (shown in Fig. [4{c)) becomes 14. 

We design a heuristic which uses a top down approach to construct the minimum 
total arrival time convergecast tree. It can be observed that while constructing 
a convergecast tree if u becomes parent of v, edge instance ef, is added in the 
convergecast tree, and w is the child of v, then we need to add edge instance e’,, 
such that tı > t. Thus the foremost journey J; from v to r in the original temporal 
graph may not be included in the constructed convergecast tree. It may need to 


choose another journey Jz such that the departure time of J2 is greater than or 
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Fig. 4. (a) A temporal graph G(V,€), (b) and (c) Two different convergecast trees on G rooted 
on node a. 


equal to the arrival time of the journey at v from u and arrival time of J2 is greater 
than the arrival time of Jı. Based on this observation we define a parameter called 
contribution of a node for a journey to the root node. 


Definition 15. Partially constructed convergecast tree: For a temporal graph 
G(V,€), if a convergecast tree T(V’, E’) is constructed by adding a non-empty set 
of nodes, S C V at a time, then T at any point of time when V Æ V’ is called a 
partially constructed convergecast tree. 


Definition 16. Reach time of node v (m): Let v be a node in a fully or partially 
constructed convergecast tree T(V’, €’). Then reach time of v, py = (mar{tileti, € 
E' and u € P(v)} + 1) where P(v) is set of parent nodes of v. If P(v) = Ø then 
by = 0. 


Definition 17. Contribution of a node v for a journey J‘? (v,7r) 
(Crests (wr) Let v be a node in temporal graph G(V,€) and t* be the arrival time 
of the foremost journey from node v to root r. Let Jt=tf(v,r) be a journey in G 
and u be the next node to v in J‘s'f(v,r). Let arrival time of the foremost jour- 
ney with departure time Hu from u to r be t*’. Then contribution of node v for 
Jéstt(u,r),C = 2tp —(t*+¢*) when u £ r and Efe 


u=r. 


= ts —t* when 


v 
J’ S (v,r) tsiti (ur) 


Figure [5]shows a partially constructed convergecast tree where r is the root node. 
The dashed lines are the edges of the temporal graph and the solid lines are edges of 
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Fig. 5. A temporal graph and a partially constructed convergecast tree 
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Algorithm 5 minTotalConvergecastTree(G(V, E), r) 
Input: A temporal graph G(V, €), lifetime T, root node r. 
Output: A convergecast tree T(V’, E’) rooted at node r, where V’ = V. 
1: computeAllJourneys(G(V, E), r) 
2: Vu € V, u.label := u.lateDep 
3: T := constructConvergecastTree(G(V, €), r) 
4: Vu E V, u.label := Hu 
5: T := constructConvergecastTree(G(V, E), r) 


the partially constructed convergecast tree. Nodes u1, ug are two parent nodes of v 
and ly is 2. The contributions of the journeys J?4(u3,r) = {((us, v), 2), (w, r), 3)} 
and J+?(us,r) = {((u3,r), 1)} are 3 and 0 respectively. 

We propose a two phase heuristic to construct a MTATCT for a given temporal 
graph G(V,€). In the first phase we label all nodes with the latest departure time 
to r. Then we add nodes to tree T(V, €’) in increasing order of labels one by one 
starting with the node with the minimum label. While adding node u, we greedily 
add edge instance e’,,, to E’ such that following conditions are satisfied: 


(1) There is a journey J‘*f(u,r) and v is the next node to u in JH (u,r). 

(2) cy =¢ or if eft., E€ E then tı >t+1. 

(3) Contribution of Jtt? (u, r) is the minimum among all available journeys from u 
to r which satisfy conditions 1 and 2. 


Node v adds u as its parent and u adds v as its child node. The first phase terminates 
when all nodes are added to T. 

In the second phase we label each node with its reach time. Then we start from 
the node with the minimum label and follow the same greedy approach to select the 
child of the node and add those nodes to the final convergecast tree. If the labels of 
two nodes are the same, we consider the node which has the higher departure time. 
The pseudocode of the algorithm is given in Algorithm [5] Each node is associated 
with some fields which are given below: 


(1) u.label: Label assigned to a node u. This label information provides the order 
in which nodes will be added to the convergecast tree. 

(2) u.lateDep: Latest departure time of journey from node u to r. Latest departure 
time for node r is T. 

(3) u.st: Timestep when journey can start from node u after addition of it to the 
convergecast tree. Initially when c, = ¢, u.st = u.lateDep. 

(4) nextNode(J‘:"'f(u,r),v): It returns the next node of v in journey J‘*"'f(u,r). 
Here v is an intermediate node in J's" (u,r). 

(5) Su: Set of foremost journeys at different timesteps from node u to r. 

(6) nodeList: List of nodes in temporal graph. 


404 S. Mandal & A. Gupta 


(a) 


Fig. 6. (a) A temporal graph, (b) Convergecast tree after first phase, (c) Final convergecast tree. 


For the temporal graph given in Fig. [6{a) with lifetime 8 and the root node 
a all non-root nodes, b,c,d,e, f,g,h are labeled with 4, 3, 2, 4, 2, 0, 7 respec- 
tively according to their latest departure time. According to the heuristic node 
g is added first to the convergecast tree T) and f gets marked as its child node 
with uf = 1. After that node f and node d are added to the tree. While adding 
d, the algorithm finds that it has two journeys J™? (d,a) = {((d,b), 0), ((b, a), 1)}, 
J?3(d,a) = {((d,a),2)} with contribution 0, 1 respectively to a. So edge instance 
ep, is added to Tı and b becomes child of d and py = 1. Then c is added with b as 
its child and u, = 4. At the end of the first phase the constructed convergecast tree 
with total arrival time 32 is shown in Fig.[6{b) where the tree edges are marked with 
the solid lines. In the first phase when node h is added to T4, though it has a journey 
J'4(h,a) = {((h, f), 1), (Cf, e), 2), ((e, a), 3)} with smaller contribution, it cannot 
be added because f is already added to Tı with departure time 1. For the second 
phase all nodes are labeled by the reach times of them. So nodes b,c, d,e, f, g, h are 
labeled with 4, 0, 0, 2, 1, 0, 0 respectively. As node h is labeled with the smaller value 
than node f, it is added before f. In the second phase when d is added to Tə the 
value of pp is 4 and the contribution of journey J°°(d,a) = {((d, b), 0), ((b, a), 4)} 
is 3 and the contribution of J>’ (d,a) = {((d,a),2)} is 1. So this time e?,, is added 
to T>. The final tree, Tz is shown in Fig. [6{c) with total arrival time 29. 


Theorem 18. Algorithm[] correctly computes a convergecast tree. 


Proof. Line[M4]of Algorithm[Jensures that for all nodes the algorithm adds exactly 
one edge instance to its child node. For any node v all the edge instances between 
v and its parents exist in earlier timesteps than the timestep in which the edge 
instance between v and its child node exists. Hence no journey from any node can 
reach the root node via its parent nodes. This proves that there is exactly one 
journey from any node to the root node. 

We have to show that there exists a journey from each non-root node to the root 
node in T by Algorithm[5] We know that any non-root node v with latest departure 
time t has a neighbour with higher latest departure time or r as its neighbour. So 
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Algorithm 6 constructConvergecastTree(G(V, €), r) 
1: Initialize: Vu € V, P(u) :=0,cy := ¢, u.st := u.lateDep, Hu := 0 
2: sort nodeList in ascending order of node.label value 
3: for i = 0 —> n do 
4: u := nodeList/i] 
5: Jést (u, r) := MinC es tp url {Jet (u, r) € Su } and ts > Hu and v.st > 
ts +1, where v := nextNode(J*=*t (u, r), u)} 


6: v := nextNode(J*=*t (u, r), u) 
T: u.st := ts 

8: ép =n 

9: 


P(v) := P(v) U {u} 
10: if wy < ts +1 then 
11: fly := ts +1 

12: end if 

13: V := V'U {u,v} 

14: 2 := E'U {ets } 

15: end for 

16: return(T’(V’, €’)) 


when node v is added by the proposed algorithm there is at least one neighbour 
through which it has a journey to the root node r. At Line D] Algorithm [6] chooses 
the journey with the minimum contribution value among all available journeys. 
Hence there is a journey to the root node r from all the nodes in T. Hence T is a 
convergecast tree. 

In the second phase, nodes are labeled with their reach times. In the tree con- 
structed after the first phase, the edge instance to child node of a node exists in later 
timestep than its reach time. Using this observation and following similar arguments 
as followed for the first phase tree, it can be proved easily that the constructed 
structure is a convergecast tree. m 


Theorem 19. Algorithm [E] computes a convergecast tree in O(T (n + m)) time for 
a given temporal graph G(V, E) with lifetime T. 


Proof. Finding the foremost journeys using Algorithm [[]takes O(T (n +m)) time. 
Sorting the nodes using labels (latest departure time in the first phase and reach 
time in the second phase) takes O(n + T) time as labels are integers and bounded 
by T. For each node there can be at most 7 number of journeys to the root node r. 
So this algorithm needs to compute contributions for at most 7 number of jour- 
neys. Thus addition of each node to the convergecast tree takes O(7) time. So the 
convergecast tree construction at each phase takes O(n +7 +7) time. Hence the 
total running time of this algorithm is O(T (n + m)). o 
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7.1. Experimental results 


We have evaluated the proposed algorithm by running it on some real life datasets 
representing different types of dynamic networks that capture the mobility patterns 
of mobile devices [7], contacts between people [16], [27], flight schedule between dif- 
ferent. airports 
and communications between different autonomous systems [22]. The dataset in 
Ref. reports the communication opportunities between 98 imotes distributed 
among the participants of INFOCOM’06. The dataset from SocioPatterns [I6], [27] 
contains contact information between students in a primary school. Our next dataset 
from includes 
flight schedule information between different airports of USA for the years 2006 to 
2008. The other dataset from [22] has 733 daily instances of communication data 
between different nodes of different autonomous systems on the inter-net. Though 
the temporal graphs formed from these datasets are non-periodic, and hence the 
changes in edge set cannot be known in advance, we considered them to see how 
the heuristic will perform on large temporal graphs. 

Each dataset contains contact information for a total of T time, where T varies 
between the data sets. For each such dataset, we have considered several values of 
p < T, and divided the total time into T = = subintervals, each of which is taken 
as a single timestep for the temporal graph to be constructed. The static graph for 
each timestep is constructed with edges added between two nodes if the nodes come 
into contact at least once in that timestep (i.e., in the subinterval of length p). This 
gives a temporal graph with lifetime 7. We have studied how different values of p 
changes connectivity between nodes. 

For the first dataset in Ref. [7], T = 342915 seconds. The values of p considered 
for this data set are 1800, 3600, 5400 and 7200 seconds. For the second dataset [I6], 
27], T = 61960 seconds. The values of p considered for this dataset are 5400, 7200, 
9000, 10800 seconds. For the flight dataset we consider each airport as node and a 
scheduled flight denotes an edge. This graph is directed in nature, so we have taken 
values of p = 1,2,3,4 months and flight between two nodes as an edge. We have 
taken T = 36 months of flight schedule as our dataset. For the dataset in [22], we 
have taken T = 50 days of communication data, with p = 1, 2,3, and4 days (T was 
taken to be 51 and 48 days for p = 3 and 4 respectively to get an integer number 
of timesteps). We have taken these values of p such that for the 1% data set the 
lifetime of the temporal graph becomes O(n) (n is the number of nodes) and for the 
204, 37d and 4t? data sets it is much less than the number of nodes. To construct 
the convergecast tree, a root node is chosen such that all non-root nodes have at 
least one journey to the root node in all four temporal graphs constructed from that 
dataset. The first nine columns in Table [] show various properties of the different 
datasets chosen. 

It can be easily seen that the total arrival time for any MTATCT of a temporal 
graph is lower bounded by the sum of arrival times of foremost journeys from every 
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Table 1. Results of running proposed algorithm on different data sets. 


Data Sets Interval Length No. of No. of Lifetime Avg no. of Avg no. of Avg hop root node Lower Obtained total 


(p) Nodes Edges (T)  instances/edge journey/node —_ count/journey bound arrival time 
1800 sec 98 4398 190 15.32 171.12 2.33 50 797 944 
eae 3600 sec 98 4398 95 10.55 85.29 2.10 50 516 937 
EOE OMG 5400 sec 98 4398 64 8.67 56.89 1.94 50 396 625 
7200 sec 98 4398 48 7.41 42.77 1.86 50 333 527 
5400 sec 242 8317 2 2.78 10.01 2.04 1838 667 768 
E TEA 7200 sec 242 8317 10 2.57 7.95 2.03 1838 556 644 
: ee 9000 sec 242 8317 7 2.18 5.64 1.96 1838 509 575 
10800 sec 242 8317 6 2.09 4.68 1.95 1838 492 530 
1 month 313 2718 36 28.51 34.54 1.47 ATL 974 2430 
EE 2 months 313.2718 18 14.53 17.06 1.44 ATL 712 1464 
lai 3 months 313—2718 12 9.83 11.22 1.42 ATL 616 1105 
4 months 313 2718 9 7.47 8.34 1.41 ATL 577 945 
1 day 3291 7124 50 41.04 48.52 2.17 3561 8331 31312 
Sas 2 days 3291 7124 25 20.82 23.47 2.17 3561 9903 48194 
BSS anne 3 days 3291 7129 17 14.26 15.59 2.15 3561 8978 34252 
4 days 3291 714 12 10.17 10.68 2.13 3561 8483 26256 


non-root node to the root node. We use this sum as a reference to compare the 
performance of the heuristic, even though it is only a loose lower bound. We have 
also studied the impact of different properties of a temporal graph on the total 
arrival time of the constructed convergecast tree by the proposed heuristic. 

The last two columns in Table [I] show the reference lower bound and the total 
arrival time of the MTATCT computed by the heuristic. It can be seen that for the 
first and second datasets, the obtained value is within two times of the lower bound 
in all cases; however, for the third dataset it is up to three times, and for the fourth 
dataset it is up to five times of the lower bound. For the first two datasets, the edge 
to node ratios are 44.79 and 34.37 respectively, for the third dataset it is 8.68, and 
for the fourth dataset it is 2.16. This shows that temporal graphs generated from the 
first two datasets are more densely connected than those generated from the third 
and fourth datasets. The arrival time of any journey from a source to a destination 
is affected by two parameters, hop count and waiting time on intermediate nodes. 
In Table [| we can see that average hop count per journey does not differ much 
between datasets. As only one journey from each node is allowed in a convergecast 
tree, any node cannot start its journey to the root node before arrival of all journeys 
from its parent nodes. Increasing the number of parents for a node v results in a 
larger difference between the arrival time of its foremost journey in the temporal 
graph and the arrival time of its foremost journey in the convergecast tree. This also 
increases the gap between the arrival times of foremost journeys in the temporal 
graph and the arrival times of foremost journeys in the convergecast tree for all 
the nodes in the subtree rooted on v. If a graph is sparse then the chance of a 
node having more number of parents increases. This explains the larger difference 
between the lower bound and the total arrival time of the computed MTATCT for 
the third and fourth datasets. 


8. Conclusion 


In this paper we have defined two types of convergecast trees for temporal graphs, 
namely bounded arrival time convergecast tree (BATCT) and minimum total arrival 
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time convergecast tree (MTATCT). A O(T (n + m)) time algorithm is proposed to 
construct BATCT. We have shown that MTATCT problem is NP-complete and 
proposed a O(T(n + m)) time heuristic for the problem. Detailed experimental 
analysis is done on four real life data sets to evaluate the proposed heuristic. 
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